Breathers in a Discrete Nonlinear Schrodinger type Model : Exact Stability Results. 
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Exact stability analysis of 1-site breathers in an NLS type model (ref. [4], see below) indicates 
destabilisation through a growth rate becoming positive as a stability border is crossed, while above 
a critical spatial decay rate (A c ) the breather is always unstable. Similar analysis on 2-site breathers 
with phase difference n (see below) indicates destabilisation through Krein collision of a pair of 
growth rates. All eigenmodes can be exactly calculated and conform to numerical estimates. 
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PACS number(s): 05.45.-a, 63.20.Pw, 63.20.Ry 

The question of dynamical stability of discrete breathers 
has been addressed in the literature and a number of use- 
ful results are now known. The stability of a monochro- 
matic breather depends on the eigenvalues of a symplec- 
tic matrix, the so called Floquet matrix. Aubry [pj in- 
troduced a novel extension to the concept of the Floquet 
matrix , thereby greatly simplifying the task of comput- 
ing the Floquet multipliers and of looking for the onset 
of instability. In particular, he established that, one-site 
breathers sufficiently close to the 'anti-integrable' limit 
are dynamically stable. However, few exact results relat- 
ing to the multipliers and the associated eigenmodes are 
known ( @, P); see references therein for a background 
to the problem) notably because the breather solutions 
around which the linearisation is to be performed are 
themselves obtained numerically, even though with high 
precision. In this context it is useful to have models that 
admit of exact breather solutions as also exact stabil- 
ity results. In Q we introduced a number of piecewise 
smooth (PWS) models including a nonlinear Schrodinger 
type one and worked out exact one-site breather solu- 
tions where one of the sites is excited to a value above a 
threshold while all other sites remain below the threshold 
( we denote these as respectively high and low). In the 
present paper we work out exact stability results for these 
breathers as also for 2-site breathers where two sites with 
an arbitrary intervening gap are high while all other sites 
are low. While the destabilisation of 1-site breathers al- 
ways occur through the growth rate of a mode becoming 
positive, the 2-site breathers get distabilised through a 
Krein collision of the growth rates (see below). 
More precisely, the PWS Schrodinger type model intro- 
duced in 101 is 



non-linear Schrodinger equation. The one-site breather 
solution centred at n = is then 
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FIG. 1. Comparison of function f(z) (a) of the usual form 
f(z) = -\z\' 2 with (b) of Eq.(2). 

The spatial decay rate A (|A| < 1) and the breather 
frequency w satisfy 

M > 2, lu = -(A + 1/A), ^ < 0, | 7 | > 1 + ±- (5) 

(note: Equations (10), (11) in Q expressing the above 
relations contained errors due to oversight, which we re- 
gret) We distinguish between breathers of type A (A > 
0,oj < —2) and type B (A < 0,uj > 2)), referring to the 
former in presenting our results below (results for type B 
are obtained analogously) . The overbar in (||) , (Q) is used 
to distinguish the breather solution from neighbouring 
perturbed solutions (see below). 

We consider these perturbations in the 'rotating frame' 
in the form 



at 

f(z) = -(l-a/\z\)Q(\z\ 



a), 



0,±1,...) (1) 

(2) 



where a(> 0), V(> 0), 7 are appropriate constants and 
denotes the Heaviside step function. We scale the 
nearest neighbour coupling strength and the threshold 
parameter to V = 1, a — 1. Fig.l compares f(z) of (0) 



with the more usual form f(z) 



of the discrete 



Ipn = (4>n + U n (t))e~ 



(6) 



Spliting u n into real and imaginary parts, u n — x n + iy„, 
one observes that the time variations of x n , y n are ob- 
tained from the Hamiltonian 

H = {{Xn+lXn + yn+lVn) + "^( x « 2 + Un 2 )} + 
n 

+ ^o 2 + yo 2 (l-J)}, 



1 



where b 



7+A- 



Eigenvalues of the matrix L deter- 



mining the evolution of the a; n 's and y n 's will be termed 
growth rates (we denote a typical growth rate by the sym- 
bol p) while a more convenient approach in the present 
context is to eleminate the y n 's or x„'s, arriving at 



X = AX, Y = AY 



(p o), 



(7) 



where we define X — ( £-2, %o> %2, ) T , 

Y = ( y_ 2 ,y-i, 2/0,2/1, Z/2, ) T , and A is a banded 

matrix with elements a mn (— oo < m,n < oo) differing 
from zero only when |m — n\ < 2 : 

a ,o = 2( 7 A-2- A 2 ), 

ao,i = ao,-i = 3A + 1/A, 

01,0 = a -i,o = 2A + 2/A - 7, 

a„,„ = -(A 2 + l/A 2 + 4), (n^O), 

«n,n+l = a n +i, n = 2(A + 1/A), (n ^ 0, -1), 

The eigenvalues (/z) of ^4 are related to growth rates 
(p) as /it = p 2 . The spectrum of eigenvalues /x and the 
corresponding eigenmodes can all be worked out exactly 
for any given parameter values A, 7 satisfying (||). There 
exists a band of eigenvalues — (ui — 2) 2 < \i < — (u> + 2) 2 
with extended eigenmodes, there being a symmetric and 
an antisymmetric eigenmode associated with each eigen- 
value. Additionally, there exists a single localised sym- 
metric eigenmode for 7 > 71(A) where 



7i(A) 



(A-i)y / ((A + ^) 2 -2(A + l)) 
(A-I) + v /((A+I) 2 -2(A+i)' 



For any given A, 71 is thus the threshold for the localised 
mode: for 7 < 71 there is no localised mode and the band 
edges are empty, while for 7 = 71 an eigenmode appears 
at the band-edge, and dissociates from the band as the 
localised mode for 7 > 71. One notes in this context 
that there always exists a trivial eigenmode (i = with 
x n = 0, y n oc (p n , which is to be distinguished from the 
eigenvectors of A, A T with fi = 0, the latter correspond- 
ing to a spurious mode in the present context. While a 
typical localised eigenvector of A is of the form 



\n\ 



n\ 



where the spatial decay rates xii Xi (each of modulus 
< 1) as also the ratio a/f3 can be obtained exactly, the 
eigenvector corresponding to [i = looks like 

x n = (a + /3|n|)A |n| . 

The growth rates for the localised mode are imaginary 
(p = ±iy[/i|) so long as 71(A) < 7 < 7(A) where the 
transition value 7 can also be obtained exactly in the 
model: 



7(A) 



(1 + 4A 2 + A 4 )(l — A 2 ) 
2A3 ' 



As 7 crosses 7 from below, the isolated eigenvalue /x 
crosses /i — and becomes positive, giving rise to a pos- 
itive growth rate (p + = this is associated with a 
negetive growth rate p_ = — yfii, making up a pair) and 
implying destabilisation of the breather. 



FIG. 2. Spectrum of eigenvalues of matrix A for (a) 
71(A) < 7 < 7(A) and (b) 7 > 7(A); cross denotes eigen- 
value for the isolated localised mode while thick black line 
denotes the band; circle corresponds to n = 0. 

Fig. 2 depicts the disposition of eigenvalues on either side 
of the destabilisation, while 
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FIG. 3. Disposition of growth rates p in the complex plane 
with (a) 7 < 7 (b) 7 > 7 

Fig. 3 indicates how the disposition of growth rates in the 
complex plane changes across the destabilisation border. 
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FIG. 4. Curves d (7 = 70(A)), C 2 (7 = 7(A)) and 
C3 (7 = 71(A)) in the A — 7 plane; the point of intersection of 
Ci and C2 corresponds to A = A c ~ 0.6948 

Fig. 4 depicts the A — 7 parameter space for the model to- 
gether with the curves C\ (7 = 7o(A) = 1 + t), C 2 (7 = 
7(A)), C3 (7 = 71(A)). The breather exists only for 
points lying above C\ and is stable only for points lying 
below C2, i.e., for (A, 7) lying in between C\ and C2 the 
model admits of a stable breather solution. The curve C3 



2 



gives the threshold for the localised mode. Interestingly, 
the model predicts the existence of a critical value (A c ) of 
the spatial decay rate corresponding to the intersection 
of Ci and C2 (A c w 0.6948) so that the breathers with 
A > A c are unstable for all values of the strength of non- 
linearity 7 (> 70(A)). All these features resemble those of 
'internal' or 'breathing' modes of DNLS breathers whose 
relevance has been emphasized in the literature ( || , || ) . 
However, the antisymmetric 'pinning' mode f|, known 
to exist for DNLS breathers, is marginal in the present 
PWS model (A = 1, fi = 0) and possesses a linear spatial 
growth. 
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FIG. 5. Numerical integration following (|l|), (Q) of an ini- 
tial profile given by (|) at t = 0; (a) A = 0.5, 7 = 4.35, 
in between Ci and C3 of Fig. 3 - profiles at values of t span- 
ning a quarter of a period (T/4) starting from 2048T; (b) 
A = 0.7 (just beyond intersection of C\ and C2 in Fig. 3), and 
7 = 2.43(> 7o (A)). 

We show in Fig. 5a the temporal evolution according to 
Eqs. ([[]), (|2|) of a profile initially (t=0) coinciding with 
(||), with (A, 7) in the region between C\ and C2 while 
Fig. 5b shows the similar evolution for a breather with 
A > A c , above C\. One notes that the breather in Fig. 5a 
performs stable oscillations while that in Fig. 5b breaks 
up within a short time, confirming our results. In this 
context see also Figs. 4,5 of 

The eigenmodes can be obtained exactly in the model. 
For instance, Fig. 6 shows the localised eigenmode for 
A = 0.5, 7 = 6 calculated from the theory as compared 
with the eigenmode obtained numerically by diagonalis- 
ing the matrix A. 



two-site breathers with various possible lengths of inter- 
vening gap A. One other feature of crucial relevance is 
the relative phase between the two sites. This relates to 
the fact that (^) admits of an arbitrary phase factor e lS 
that becomes relevant in 2-site or multisite breather so- 
lutions. Thus, considering a breather with two sites high 
(i.e., with \ip n \ > 1) and all the other sites low (\ip n \ < 1) 
with the high sites having opposite phases (6 — tt) one 
obtains a monochromatic breather solution of the form 
(^) with the breather profile now given by 

J. - b r (A»-A w -") 



1-X N 
= b\~ n 
= -b\ n - N 



(0 < n < N) 

(n < 0) 
(n > AO, 



(8) 



where b 

l+A 



7(1-*") 



7 (i_aw)-(a- 1 -a) an< ^ ^' 7 are t° satisfy 
7 > jzttn ■ Here the two high sites have been chosen 
to be at n = and n = N for any given separation 
N. It is now easy to construct the stability matrix A 
of Eq.(0) and to obtain the spectrum of eigenvalues to- 
gether with the exact eigenmodes. There is once again a 
band extending from — (uj — 2) 2 to — (ui + 2) 2 (recall that 
we are considering only type A breathers, with lu < —2) 
containing extended eigenmodes, together with a number 
of isolated eigenvalues corresponding to localised modes. 
The latter appear above a certain threshold value of 7 
(for any given A), and are now typically four in number 
of which one is the trivial eigenvalue /i = while the oth- 
ers lie in the range — (u> + 2) 2 < fi < 0. Ordering these as 
Mi < M2 < M3 < M4 — 0j t ne eigenmodes corresponding 
to fix and /X4 are antisymmetric (with respect to the sites 
n = 0, A) while the other two are symmetric. Fig. 7 shows 
the four eigenmodes for A = 8, 7 = 4.9443, A = 0.5. For 
a given A, A, as 7 is made to increase from the thresh- 
old value, a boundary 7// (A) is reached when /12 and /X3 
become coincident and then become complex. In terms 
of the growth rates p (p — i-y/Ji) there occurs a Krein 
collision of a quartet of growth rates on the imaginary 
axis and the subsequent spliting of the growth rates out 
of the complex plane, as depicted schematically in Fig. 8. 



y=6, h=0.5, n=-0.0495 




FIG. 6. Eigenmodes obtained (a) from exact analysis and 
(b) by diagonalising the matrix A 

As we mentioned in H, the PWS model ([!]), (^) is a 
veritable little 'laboratory' yielding exact breather solu- 
tions of a wide variety. For instance, one can construct 



gap between peaks = 7 

^4.9443, fc=0.5 

H.=-0.2501 




u=-0.1154 



U=0 



site number (n) 

FIG. 7. Symmetric and antisymmetric localised eigen- 
modes for a 2-site breather. 
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In a model with sufficiently smooth f(z) (ref. Eq.||) this 
leads to the Hamiltonian Hopf bifurcation ( [Q , see refer- 
ences therein; the Krein collision has been discussed in 
in terms of the Floquet multipliers related to the breather 
solution - the latter are related to our eigenvalues (n) and 
the corresponding growth rates (p) in a simple manner; 
there exists a large body of literature relating to Krein 
collision in symplectic mappings - see, e.g., ||). In the 
present context, the collision signifies a loss of stability. 
For instance, with N — 8, A = 0.5, the model predicts a 
loss of stability of the two-site breather through a Krein 
collision at j H ~ 4.9444. 

Fig. 9 depicts the time evolution of an initial profile cho- 
sen in accordance with (^) with 7 on either side of the 
collision: one notes that the breather indeed remains sta- 
ble for 7 < jh and breaks up quickly for 7 > 7^. 
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high sites are then almost unaffected by each other, each 
essentially behaving like a one-site breather (3), (4). 
Interestingly, a 2-site breather with the two high sites 
in phase with each other can also be constructed in the 
model, but turns out to be always unstable. 
The details of the calculations together with numerical 
and other results on breather stability will be presented 
elsewhere. 
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FIG. 8. Sequence of changes in the disposition of growth 
rates in the complex plane; thick black lines denote the band 
while the dots and circles correspond to localised modes; (a) 
before Krein collision^ < "/h; A fixed), (b) Krein collision 
7 = jh, (c) a quartet of complex growth rates resulting from 
the Krein collision, and a pair of imaginary growth rates ap- 
proaching collision at p = 0, (d) Krein collision 'in the re- 
verse', (e) all four growth rates of localised modes real. 



FIG. 9. Numerical integration following (1), (|2j) of an ini- 
tial two-site breather profile given by (^): (a) 7 < 7h (A), (b) 
7 > Jh{\). 



As mentioned above, in Fig. 8 we depict the sequence of 
changes in the disposition of the growth rates p as 7 is 
made to increase through 7^ for fixed A. One observes 
that at some value 7(> jh) of 7, the root /i becomes posi- 
tive (this is similar to the crossing of a negative eigenvalue 
through the value zero at 7 = 7(A)) with \i2, M3 remain- 
ing complex; subsequently at some value 7 = j' H the lat- 
ter collide at a point on the positive real axis and move off 
in opposite directions - Krein collision 'in the reverse'. In 
the present context of stability of the breathers, however, 
this reverse collision does not have a direct relevance since 
the breather, having become unstable at 7 = jh contin- 
ues to remain unstable even at 7 = 7 and 7 = 7^-. For 
large N, one finds that 7//, 7 and j' H are all contained 
in a small interval around 7 as expected, since the two 
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